Fermionization in an expanding ID gas of hard-core bosons 
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We show by means of an exact numerical approach that the momentum distribution of a free 
expanding gas of hard-core bosons on a one-dimensional lattice approaches to the one of noninter- 
acting fermions, acquiring a Fermi edge. Yet there is a power-law decay of the one-particle density 
matrix p x ~ 1/y/x, as usual for hard-core bosons in the ground state, which accounts for a large oc- 
cupation of the lowest natural orbitals for all expansion times. The fermionization of the momentum 
distribution function, which is not observed in equilibrium, is analyzed in detail. 
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At very low temperatures and densities a one- 
dimensional (ID) gas of bosons is expected to behave 
as a gas of impenetrable particles known as hard-core 
bosons (HCB) Two recent experiments successfully 
achieved the requiredparameter regime and made HCB 
a physical reality @, ■ In contrast to bosons in higher 
dimensions, ID HCB share many properties with nonin- 
teracting spinless fermions to which they can be mapped 
. Thermodynamic properties like the total energy, and 
microscopic properties like density profiles are identical 
in both systems. On the contrary, quantities like the mo- 
mentum distribution function (njt) [f| and the so-called 
natural orbitals (NO) || are very different for HCB and 
spinless fermions. This is due to the different behavior 
of the off-diagonal elements of the one-particle density 
matrix (OPDM) in both systems (see e.g. Ref. Q). 

An important point for the comparison between ex- 
perimental results and theory in Refs. 0, has been 
the awareness of the effects of the trapping potential in 
the properties of the HCB gas. On a harmonic trap it 
has been found that the power-law decay of the OPDM 
Pij ~ \xi — Xj\~ 1 ^ 2 , known from homogeneous systems 0, 
is renormalized by a factor that depends on the density at 
points i and j. This factor is proportional to [riirij} 1 ^ 4 for 
continuous systems 0, and to [n,(l — ni)nj(l — "■j)] 1 ^ 4 
for HCB on a lattice Furthermore, the power-law 
decay of the OPDM has been found to be universal in- 
dependently of the power of the confining potential 0- 
Even for systems out of equilibrium, that start their evo- 
lution from a totally uncorrelated state, the power law 
above develops dynamically and leads to the emergence 
of quasi-condensates at finite momentum (lpj . 

In this work we show that during the free expansion 
of ID HCB in a lattice a further degree of fermionization 
takes place, i.e., rik of the HCB becomes for long expan- 
sion times equal to the one of nonintcracting fermions, 
displaying a Fermi edge. This feature, absent in equilib- 
rium, can be easily confirmed experimentally with a setup 
like the one in Ref. [j| where the expansion of the ID gas 
was studied after removing the axial confinement. Other 
quantities like the OPDM and the NO still evidence the 
strongly interacting character of the HCB system. 



We obtain the exact dynamics of expanding clouds of 
HCB on ID lattices on the basis of the Jordan- Wigncr 
transformation, which maps the HCB Hamiltonian into 
the one of spinless fermions . Apart from being exact, 
our method allows to consider relatively large number of 
particles and system sizes, which can be even larger than 
the ones in the present experimental setups. Results for 
continuous systems can be extrapolated from very low 
densities in the lattice At time t = 0wc switch off 
the trapping potential, and start the time evolution from 
the ground state of the HCB Hamiltonian 
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which has the additional on-site constraints b\ 2 = bf = 0, 
{&i,&t} = 1. InEq. (J), b\ and h represent the HCB cre- 
ation and annihilation operators, respectively, rij = b\b, t 
the particle number operator, t the hopping parameter, 
and Vi the curvature of the harmonic confining potential. 

In order to characterize the initial state, we use the 
characteristic density p = N^a/C, 0, [^}, which for 
trapped systems plays a similar role than the density 
in periodic systems. Nt is the number of HCB, £ = 
0^2 A) is a length scale of the trap in the presence 
of the lattice, and a is the lattice constant. The quasi- 
momentum distribution function is also normalized 
by the length scale ( as rik — (a/Q J2ij e ~ ~ Piji w hh 
Pij = (b\bj). In addition, we define the Fermi momen- 
tum hp associated to the fermions as ep = — It cos(kpa), 
where ep is the energy of the last occupied fermionic 
single-particle state in the trap at t = 0. (Although 
in the trap is continuous at hp, it is possible to see 
that rife approaches zero even faster than exponentially 
for k > kp.) 

In Fig. n we show for an expanding gas of HCB 
at four different times, and compare it with the one of 
nonintcracting fermions (which does not change during 
the expansion). Several issues are evident: i) Shortly 
after switching off the trapping potential, the peak at 
n fc=o disappears, ii) For k < kp a redistribution of the 
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population of k states takes place, such that starting from 
k ~ 0, rife of HCB matches in time the one for fermions. 
iii) k states with k > kp become less populated and an 
edge develops at kp. The overall process leads to an 
for the HCB that is equal to the one of the fermions. 




FIG. 1: (color online), for 100 HCB expanding from an 
initial state with p — 0.51, compared to the one for the cor- 
responding fermions. Times (r) are given in units of h/t, and 
fci? denotes de Fermi momentum as defined in the text. 

After observing the fermionization of rik for HCB, one 
could naively expect that something similar could be hap- 
pening with the NO occupations. The NO (</>^) are the 
eigenfunctions of the OPDM J2j Pijtf] = \4>i El> i- e -> 
they are effective single-particle states with occupations 
A I( . For noninteracting fermions the NO are the eigen- 
functions of the Hamiltonian and their occupation is one. 
We find that in contrast to rife the NO occupations do not 
fermionize, as can be seen in Fig. [21 for the same param- 
eters of Fig. n 
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FIG. 2: (color online). Occupation of the NO vs r\ for the 
same initial trap parameters and times of Fig. The thin 
dashed line corresponds to a power law r]~ 4 , which is known 
from equilibrium systems at low densities [j| . 

There are two features of the NO occupations that we 
find worth noticing. The first one, which is difficult to 
distinguish in Fig. [2 is that the lowest natural orbitals 
slightly increase their occupations during the expansion 
of the gas. This can be intuitively understood as an in- 
crease of the "coherence" of the system due to an increase 
of the system size, which delocalize the HCB over more 
lattice sites. Something similar occurs in the ground state 



occupations of the lowest natural orbitals when for the 
same number of particles the curvature of the trap is de- 
creased 0- However, in equilibrium systems, n k=0 also 
increases with the increase of Ao- The apparent contra- 
diction between the decrease of nk=o and the increase of 
A in the expanding gas can be resolved observing the 
Fourier transform of the lowest NO at r = and t > 0. 
As it can be seen in Fig. initially \</)®\ has a peak at 
k = showing that quasi-condensation occurs around 
k = 0, and this is reflected in rife. For r > the lowest 
NO becomes an extended object in fc-space so that the 
HCB forming the quasi-condensate have many different 
momenta, basically as many as in Fig.^ Hence, there 
is no contradiction between the observed behavior of the 
NO occupations and n^, although the last one is clearly 
different to the one of systems in equilibrium. 
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FIG. 3: (color online). Fourier transform of the lowest NO, 
and rik in a logarithmic scale (inset), for the same initial trap 
parameters and times of Figs. Q and |21 &f (vertical dashed 
line) denotes de Fermi momentum. 

The second feature that is worth noticing, sets in only 
when the density of the expanding HCB becomes very 
low. Then the universal power-law decay of the NO occu- 
pations for larffe-77 (A r) ~ AM b ri~ A )., known from equilib- 



rium systems also appears in noncquilibrium (Fig.[2*)l. 
In Ref. 0, the prefactor of the power-law Ajj b was found 
to depend only on Nb independently of the confining po- 
tential. We find out of equilibrium that Ajv 6 is exactly 
the same than in the ground state case ||. For Ni, = 100 
we have plotted A r; = An^ 4 in Fig. [2 The long tail 
of the momentum distribution function ~ |&| -4 , that 
in equilibrium appears together with the A,, ~ r\~ 4 0, is 
in general not present in noncquilibrium since for large 
r, rift for HCB starts to behave like the one of fermions 
(inset in Fig.^J. 

Considering the previous results for the behavior of the 
NO occupations, which is similar to the one known in 
equilibrium 0, one expects that also the OPDM should 
behave similarly. Since in nonequilibrium pij = \pij\e ' ij 
is in general a complex object, in order to compare with 
equilibrium systems we first study its modulus. Results 
for the same systems of Figs. ^Glare shown in Fig.0{a). 
Figure Ufa) shows that |py (r)| have exactly the same 
form than pij in equilibrium systems || . For large values 
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of \xi — xj\ a power-law decay \pij\ ~ \xi — Xj\~ x l 2 can 
be observed for all times, and the prefactor of the power 
law decreases with the reduction of the local densities in 
the system. Hence, the slow decay of the one-particle 
correlations is the one accounting for the large (~ \/Nb) 
occupation of the lowest NO. On the other hand, Figs. 
EHb)-(d) show that the phase of the OPDM (%) starts 



to increasingly oscillate at large distances. In particular 
Fig. 0Jb) shows that after a very short time, when the 
modulus of the OPDM have almost not changed, 6*^ have 
started to oscillate for \xi — Xj\ ^> a producing a fast 
destruction of the zero momentum peak in n b k , as shown 
in Fig.Hl 
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FIG. 4: (color online). Modulus of the OPDM (a) and its phase (b)-(d) for the same initial trap parameters and times of Figs. 
0E3 Both quantities have been evaluated with respect to the middle of the expanding cloud of HCB. Thin continuous lines in 
(a) correspond to power laws \xi — Xj 
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In order to gain further insight into the fermionization 
described above, we observe that for noninteracting par- 
ticles, given the initial one-particle density matrix p a &, 
the density profile at time r can be evaluated as [llT| 



Tli(T) 



a b 



G i, a ( T ) G i,b(T)p ab , 



(2) 



where the free one-particle propagator can be written in 



the form G^fr) = J^k 



-ir / h[ek—hk(xi — Xj) / r] 



e k being 



the dispersion relation. Fourier transforming the initial 
OPDM into Qk!,k 2 an d performing straightforward inte- 
grations, Eq. can be rewritten as 



n,-(r) 



E 



ir/h[e kl -hkix/r] -ir/h[e k2 -hk 2 x/r] 



re 



k—ko (xi / t) 



~ Qk\—ko{xi I t) ^k2—ko{x.i I r) i 



(3) 



where in the last step we have assumed r to be very large 
and performed the summations using the saddle point 
approximation, [kofei/r) is determined by the expres- 
sion e' ko = tiXi/r, where the prima means fc-derivative.] 
Equation shows that for r — > oo the density profile 
is only determined by a rescaling of the diagonal part of 
0k!,k 2 j which is nothing but n^. In the continuum limit 
ft = h 2 k 2 /2m and n{x,r) ~ (m/r)n(fc = mx/Hr) [n(k) 
is the momentum distribution function] , which is a known 
result that can be also obtained by other means pd| . 
Since the arguments above are valid for both fermions 
and bosons, the fact that of HCB converges to that 
of fermions after expansion could be explained at this 



point if after a certain time, the expanding HCB could 
be considered as non-interacting. This would mean that 
rik for HCB would be determined by a rescaling of the 
density profile that, on the other hand, both in the con- 
tinuum |l2l | and on the lattice 0, ^} is the same as that 
of fermions due to the mapping connecting both. 

However, a non-interacting treatment of the HCB ex- 
pansion is invalidated by Fig. 4(a). There it is shown that 
at all times, even after fermionization, the density ma- 
trix decays as 1/y/x, a power that corresponds to HCB, 
and hence, the system corresponds to strongly interact- 
ing particles. Therefore, the expansion out of equilibrium 
leads to a new kind of bosonic state, with a Fermi edge in 
the momentum distribution function but still the effec- 
tive one-particle states, as given by the natural orbitals, 
exhibit a high occupation as expected for bosons. 

In what follows we study how exactly the fermioniza- 
tion of the momentum distribution for bosons n\ oc- 
curs in time, and how it depends on experimental pa- 
rameters like number of particles and characteristic den- 
sities. For that we analyze the relative area between 
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■ \)/Q2k n k)' This is shown in Fig. Eta) 



for p = 0.51 and different fillings of the trap. Fig. Eta) 
shows that changes in n b k occur fast in terms of the char- 
acteristic time of the system, which is given by H/t. In 
addition, if one chooses a criterion like <5 = 0.05 to state 
that n\ has fermionized, it is possible to see in the in- 
set of Fig. OJa) that for a given characteristic density the 
fermionization time (tf) grows linearly with the number 
of HCB in the trap. 
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Another question that is important to answer is the 
consequence of changing p in the ground state of the trap. 
In order to compare systems with different p, i.e., differ- 
ent rik, we display in Fig. |SJb) the ratio R between the 
size of the cloud once S = 0.05 and its initial size. Fig. 
E{b) shows that with decreasing p the ratio R reduces 
up to ~ 2.5, and that for p > 0.5 it increases very fast. 
For low p, such that the averaged interparticle distance 
is much larger than the lattice spacing, the initial lattice 
gas is equivalent to the one in continuous systems. This 
means that a fermionized n k will be more easily observed 
in continuous systems Q than in the lattice Q- (In the 
continuous case, the asymptotic fermionization of n b k was 
obtained previously in Ref. [13 •) ^ n addition, the inset 
in Fig.Efb) shows that the ratio R remains basically con- 
stant for a given characteristic density when the number 
of particles in the trap in changed. 




FIG. 5: (color online). Fermionization of n\ during the ex- 
pansion of the gas. (a) Decrease of S (see text) as a function 
of time for p = 0.51; Nb = 51 (dashed line), Nb = 100 (con- 
tinuous line), Nt = 150 (dotted line), and Nb = 200 (dashed- 
dotted line). The inset shows the fermionization time tf (see 
text) vs Nb, for p = 0.51. (b) Ratio R between the size of gas 
for 8 = 0.05 and its original size vs p, Nb = 100. The inset 
shows R vs N b for p = 0.25 (x) and p = 0.51 (V). 
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FIG. 6: (color online). Fermionization of n\ during the ex- 
pansion of the gas. For p — 3.0 and 4.0 a Mott insulator is 
formed in the middle of the trap. The behavior for a pure Fock 
state (Fock) with one particle per lattice site is also shown. 

that still the dynamically generated quasi-condensates, 
which have k ~ ±tt/2, could be used as atom lasers [Toj . 

In summary, we have shown that during the expansion 
of a ID gas of HCB the momentum distribution function 
becomes equal to the one of the equivalent noninteracting 
fermions. This is an effect that can be seen experimen- 
tally in systems with Q and without Q an optical lattice 
along the ID axes. On the other hand quantities like the 
NO and the OPDM still display the known behavior in 
equilibrium systems. In this way starting from a strongly 
interacting ID Bose gas one can realize a very unconven- 
tional system of bosons displaying a Fermi edge on its 
momentum distribution function. 
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Increasing the characteristic density of the initial sys- 
tem beyond the values in Fig. EJb) one starts observ- 
ing a behavior of S which is different to the one seen in 
Fig. Efa). The reason is that particles become more lo- 
calized in the middle of the trap, and eventually after 
p ~ 2.6 — 2.7 a Mott insulator appears in the system. 
This localization effect also generates an n\ which ap- 
proaches to the one of the fermions in the initial state. 
[In the limit of all lattice sites with occupation one [ic| . 
n k( T = 0) = n k .] When such systems are released from 
the trap, quasi-long range correlations start to develop 
between initially uncorrelated particles and they lead to 
the formation of traveling quasi-condensates [lOj. This 
generates an n k (r) at short times that may be more dif- 
ferent to n k than n|(r = 0), as it can be seen in Fig.0 
However, after long times one can see that a fermioniza- 
tion of n b k starts to occurs as before for smaller p. One 
should notice that as shown in Fig. El the time scales for 
the fermionization process for large p are very long, and 
always start affecting the low-momenta region first, so 
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